Computation of inflationary cosmological perturbations in chaotic inflationary 
scenarios using the phase-integral method 
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Abstract 

The phase-integral approximation devised by Froman and Froman, is used for computing cosmological perturbations in the 
quadratic chaotic inflationary model. The phase-integral formulas for the scalar and tensor power spectra are explicitly obtained 
up to fifth order of the phase-integral approximation. We show that, the phase integral gives a very good approximation for 
the shape of the power spectra associated with scalar and tensor perturbations as well as the spectral indices. We find that 
the accuracy of the phase-integral approximation compares favorably with the numerical results and those obtained using the 
slow-roll and uniform approximation methods. 
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I. INTRODUCTION 



The study of the spectrum of anisotropics of the cosmic microwave background radiation and inhomogeneities in 
the large scale structure of the universe provides key elements in the study of the early universe. The results reported 
by WMAP favor inflation [J, i, i] over other cosmological scenarios. The present and future observational data will 
permit us to validate and discriminate among different inflationary models. According to WMAP5 data the power-law 
inflation, the hybrid inflation and quartic chaotic inflationary models A(/>^ are ruled out, while the quadratic chaotic 
inflationary model m^t/)^ agreed with the observational data [H, 01 . 

In order to compare with observations, we should be able to obtain very accurate results for the predicted power 
spectrum of primordial perturbations for a variety of inflationary scenarios. In general, most of the inflationary 
models are not exactly solvable and approximate or numerical methods are mandatory in the computation of the 
scalar and power spectra. Traditionally, the method of approximation applied in inflationary cosmology is the slow- 
roll approximation which produces reliable results in inflationary models with smooth potentials, but cannot be 
improved on a simple way beyond the leading order. Recently, some authors have applied alternative approximations, 
such as the WKB method with the Langer modiflcation [1,13, HI, the Green function method [9], and the improved 
WKB method [l3JIl|l- 

Habib et al |l2l ll3l. 1^ have successfully applied the uniform approximation method in the calculation of the scalar 
and tensor power spectra and the corresponding spectral indices for the quadratic and quartic chaotic inflationary 
models, showing that the the uniform approximation gives more accurate results than the slow-roll approximation. 
Casadio et al [15] have applied the method of comparison equation to study cosmological perturbations during 
inflation. The comparison method is based on the uniform approximation proposed by Dingle fl9| and Miller [l7j and 
thoroughly discussed by Berry and Mount p^ . 

Recently [l^, [l^] , the phase-integral approximation [111, [H, has been applied to the calculation of the power 
spectra and spectral indices in the power-law inflationary model, showing that the phase-integral method gives results 
which are comparable or better than those obtained using slow-roll or the uniform approximation. It is purpose of 
this paper to compute approximate solutions for the scalar and tensor power spectra and their corresponding spectral 
indices for some chaotic inflationary models with the help of the phase-integral approximation method. We show that 
for the quadratic chaotic inflationary model the fifth-order phase integral approximation gives more accurate results 
than those obtained using the WKB or uniform approximation methods. For the quartic chaotic inflationary model 
we also obtain very accurate results but we do not report them since this model is ruled out by the observational data 

Q- . . 

The article is structured as follows: In Sec. II we apply the phase-integral approximation to the chaotic quadratic 
inflationary model, we numerically solve the equation governing the scalar and tensor perturbations and compare the 
results for the power-spectra obtained using the phase-integral approach with those computed with the slow-roll and 
uniform-approximation methods. In Sec. Ill we summarize our results. 



II. PHASE-INTEGRAL APPROXIMATION FOR THE POWER SPECTRUM IN THE CHAOTIC IN- 
FLATION 



In this section we discuss the application of the phase-integral approximation to the computation of the power 
spectrum in the chaotic inflationary ^m'^cj)'^ model. We apply the phase-integral approximation in the study of the 
evolution of the mode k equations for the scalar and tensor perturbations in order to compute the scalar and tensor 
power spectra. Using the flfth-order phase integral approximation we compute the scalar and tens or p ower spectra 
and their corresponding spectral indices. A detailed description of the method is given in reference |l9l. [20l|. 



A. The model 

The chaotic inflationary model was introduced by Linde [2^ [25| , he proposed that the preinflationary universe was 
chaotic which means that the fields would take different values in different points of the space following a random 
pattern and inflation will occur in virtually any universe that begins in a chaotic, high energy state and has a scalar 
field with unbounded potential energy. The simplest form of the inflaton potential V{4)) in a chaotic model is given 
by the quadratic potential 

VicP) = ^m^(b\ (1) 

giving as a result a free scalar field with mass m. 
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B. Equations of motion 



In an inflationary universe driven by a scalar field, the equations of motion for the inflaton 
parameter H are given by 



and the Hubble 



3H(f>- 



dV{(t)) 



(2) 



(3) 



where the dots indicate derivatives with respect to physical time t. In the quadratic chaotic inflationary model the Eqs. 
([2]) and ([3]) are not exactly solvable in closed form; they can be solved numerically or using the slow-roll approximation. 
In the slow-roll approximation [2^ we consider that the scalar held V{(j)) varies very slowly ^ y{4>)- Using this 
approximation, we obtain that Eq. ^ and Eq.© for the quadratic chaotic inflationary model reduce to [131 



3ff0 
H 



mcj) 



(4) 
(5) 



Using Eq. ^ and Eq. ([5]) we obtain that, in the slow-roll approximation, the expansion factor a{t) and the inflaton 
field (t){t) are 



(6) 

ai exp [ ^ ^ ; t —t^ ) , (7) 
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where 4>i is a constant of integration corresponding to the initial value of the inflaton. Eq. ([7]) shows that the Universe 
expands exponentially during inflation The slow-roll parameter is Q 



ei 



H 



(8) 



The inflationary epoch finishes when ei = 1, that is (j)f = V2Mp\, when the scalar field starts to oscillate. The mass 
m of the inflation can be fixed using the amplitude of the density fluctuations detected by WMAP. In order to fit m 
to the observational data, we demand that m ~ 10~^ Q. 

The equations of motion Q and ([3]) are numerically integrated in the physical time t. We solve the system of 
coupled differential equations ^ and ^ with the help of the sixth-order Runge-Kutta method [1^, which can be 
written as 



V6M] 



PI 



2V{cl>) 



-0, 



(9) 



1/2 



= 0. 



We choose the initial value of the inflaton as 4^ = 15.4 Mp\. Since the evolution of the inflation cf) is governed by 
a second a second-order differential equation, we need to fix the initial value for the velocity of the scalar field 4'i, 
which can be obtained using the slow- roll approximation The initial value for Ui is chosen as = 1, the mass of 
the infiaton = 1.89 x 10^^^ ^^pi- The initial condition has been selected in order to guarantee enough inflation. 
In order to find the number of e-folds, we rewrite the system of differential equations ^ as 



+ SNcp H — — = 0, 

0(j) 

N = V&Mpi hv {(/)) + (p- 



(10) 



1/2 
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We find that the inflation finishes a.t tj = 1.30 x 10' Mp^ = 3.51 x 10 ''°s, a result that corresponds to 59.84 e-folds 
before the scalar field starts to oscillate. Fig. [T] shows the evolution of the scalar field (p. 
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FIG. 1: Evolution of the potential (f) for the chaotic inflationary model im^(^^ 



In order to apply the phase-integral approximation to higher orders it is necessary to calculate the integrals uj{z) 
in the complex plane, therefore, it is of help to have an analytic expression for a{t) and (l)(t) which can be obtained 
after fitting the numerical data. We find that a(t) and (t){t), take the form 



a{t) 

m 



Gi exp (at — f3t^^ 
1 - crt, 



where a = 8.6551 x lO^^Mpi, /? = 3.1380 x IQ-^^ M^^, 7 = 15.3992 Mpi and a = 1.1201 x 10 
of the expressions pT|) and IT^ we obtain that zs is given by: 



(11) 
(12) 

With the help 



^s(i) 



{a - 2(3t) 



exp [at - (3t^ 



(13) 



Fig. |2(a)| and Fig. |2(b)| compare the fitting with the numerical result and the slow-roll approximations for a{t) 



and 4>{t). The fitting is valid up to t = 5.00 x lO^Afpj^. The inset is an enlargement of the figure. Fig. 3(a) and 



Fig. |3(b)| show the ratio of the fit and slow-roll to exact solution for a{t) and 0(t) and observe that the fitting better 
approximates the numerical result than the slow-roll approximation, therefore the expressions for the scalar and tensor 
perturbations will be constructed using the fitting a(t) and zs{t) given by expressions (|lip and (|13p . respectively. If 
we use Osr and 0sr in order to calculate the power spectrum, the expression we obtain does not approach the exact 
result. 



C. Equation for the perturbations 



Since the expansion factor a and the field cf) exhibit a simpler form in the physical time t than in the conformal time 
T], we proceed to write the equations for the scalar and tensor perturbations in the variable t. The relation between t 
and 77 is given via the equation dt = a drj. In this case, the equation for the perturbations can be written as 
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FIG. 2: (a) Evolution of the scale factor a and (b) evolution of the inflaton </> for the chaotic inflationary model ^m^4>^. Solid 
line: numerical solution; dashed line: fitting; dotted line: slow-roll approximation. The inset is an enlargement of the figure. 
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FIG. 3: (a) Ratio of the fit and slow-roll to the exact solution for a{t): ai/a and (b) ratio of the fit and slow-roll to the exact 
solution for (f>{t): 0i/0, in the chaotic inflationary model ^m^(j)^. Dashed line: fitting; dotted line: slow-roll approximation. 
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a . 

Uk + -Uk + 
a 



(az's + az's) a 



Uk 



v'k + -Vk + Tfc^ - (d^ + aa)] Vk = 0. 



(14) 
(15) 



In order to apply the phase-integral approximation, we eliminate the terms iik and Vk in Eq- l|14p and Eq. (fT5 



We make the change of variables Uk{t) 
equations: 



=^ and Vk{t) = , obtaining that Uk and 14 satisfy the differential 



with 



Uk + Rs{k,t)Uk = 0, 
Vk + RTik,t)Vk = 0, 



Rsik,t) = ^ 



(dfs + az's) a 



zs 



^T(fc,t) = ^[k^-{d^+ati)]+^{a'-2c 



(16) 
(17) 



(18) 
(19) 



where U{k) satisfies the asymptotic conditions 



Uk 
Uk 



Ak\/ a{t)zs{t), kt ^ oo 



a(i) 
2fc 



exp [— ifc?7(t)], fci— >0, 



(20) 
(21) 



the asymptotic conditions (|20p and pTjl also hold for Vk- 

We now proceed to write the explicit equations for quadratic chaotic inflation. From Eq. and Eq. p^ . with 
Eq. ([TT|) and Eq. ((12]) we obtain 



fc2 

Rs{k,t) = — exp [-2t(a - /3t)] 



32/32 + 9 (a - 2(ity 



4:{a-2(3tY 



3f3, 



1,2 n 

RT(k,t) = -^cxp\-2t(a- pt)]- - (a-2(3tf + 
af 4 



In order to apply the asymptotic condition (|2ip. we use of the relation between rj and t, which is given by: 



V 



2a,^ 



TT I —a 

exp 



4/3 



Erfi 



-a + 2fit 



Erfi 



a - 2/3to 



(22) 
(23) 

(24) 



where Erfi(z) is the imaginary error function 29]. Since the conformal time rj is defined up to an integration constant, 
the lower limit U of the integral 



dr] 



' dt_ 



(25) 



is chosen in order to make 77 = at the end of the inflationary epoch, i.e, ti = 10^ -^^pi"'^- The dependence of 77 on t is 
shown in Fig. [H We can observe that as 



-kri^O =^ kt^ 00, (26) 
-kr]-* 00 ^ fci -> 0. (27) 

Eq. (fT6)) and Eq. p?|) . where Rs{k,t) and RT{k,t) are given by Eq. ((22|) and Eq. ([23|) . do not possess exact analytic 
solution. In order to solve the differential equations governing the scalar and tensor perturbations in the physical 
time t, we use the fifth-order phase integral approximation and compare this results with the slow-roll and uniform 
approximation. 
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FIG. 4: Behavior of t; as a function of the physical time t for the chaotic inflationary model 



D. Phase-integral approximation 



In order to solve Eq. and Eq. pT]) with the help of the phase-integral approximation, wc choose the following 
base functions Q for the scalar and tensor perturbations 



QUk,t) 
QT{k,t) 



RT{k,t), 



(28) 
(29) 



where Rs{k,t) and RT{k,t) are given by Eq. (|22p and respectively. Using this selection, the phase-integral 
approximation is valid as /ct — + cxd, limit where we should impose the condition (j20[) . where the validity condition 
/i <g; 1 holds. The selection, given in Eq. I|28p . makes the first order phase- integral approximation coincide with the 
WKB solution. The bases functions Qsik,t) and QTik,t) possess turning points tj-d = ts = 1.38096 x 10^ Mpj^ and 
tret = tt = t = 1.38196 X lO^Mpj"'^, respectively for the mode k = 1.369 /i Mpc^^. The turning point represents the 
horizon. There are two ranges where to define the solution. To the left of the turning point < t < tret we have 
the classically permitted region Qgr^ikji) > and to the right of the turning point t > t^d corresponding to the 
classically forbidden region Qg j,{k, t) < 0, such as it is shown in Figs 5(a) and Fig. 



6(a) 



The mode k equations for the scalar an tensor perturbations (jl6[) and (|17p in the phase-integral approximation has 
two solutions: For < t < t^ct 




%'\k,t) 


COS 


|t^s(fc,t)| - 


TT' 

4. 


%"\k,t) 


cos 


\uJs{k,t)\ + 


TT" 

4. 


q^^^^{k,t) 


COS 


\cJT{k,t)\- 


TT" 

4. 


q^^^^{k,t) 


COS 


\uJT{zk, t\ 4 


TT" 

' 4. 



(30) 



(31) 



and for t > U, 
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FIG. 5: (a) Behavior of the function Q\{k,i). (b) Contour of integration rT-g(t) for < i < ts- (c) Contour of integration 
rrs(t) for t>Ts- The dashed line indicates the part of the path on the second Riemann sheet. 



+ 



Cl 



di 
2i/a(t) 



% ^'^{k,t) exp[- \ws{k,t)\] 
{k,t) exp[|a;s(A;,i)|] , 

exp[- |a;T(A;,i)|] 

<lT^''^{k,t) exp [\ujT:{k,t)\] . 



(32) 



(33) 



Using the phase- integral approximation up to fifth order (2A'' + 1 = 5^ N = 2), we have that qs{k, t) and qT{k, t) 
can be expanded in the form 



9 




FIG. 6: (a) Behavior of Q^{k,t). (b) Contour of integration rTT^{t) for < t < tt. (c) Contour of integration rTT^(t) for t > tt- 
The dashed lined indicates the part of the path on the second Riemann sheet. 



qsik, t)^J2 ^2ns (k, t)Qs{k, t) = [Fos {k, t) + ¥2^ (fc, t) + Ya^ (fc, t)] Qs{k, t), 

n=0 
2 

qT{k,t) = ^r2„^(fc,t)gT(fc,t) = [Yo^{k,t) + Y2^{k,t) + YA^{k,t)] QT{k,t). 



(34) 
(35) 



n=0 



In order to compute qsik, t) and qT{k, t), we compute [k, t), Kjg (fc, t), (^j t), Y^^ (k, t) and the required functions 
eQf^{k,t), e2s{k,t), eo^{k,t) and e2T^{k,t). The expressions and (|35p give a fifth-order approximation for qs{k,t) 
and qT{kjt). In ord er to compute a;s(A:, <) and tt'T(fc, t) we make a contour integration following the path indicated in 
Figs. [5(b)]- (c) [6(b)} (c) 
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ws(A;,t) = uJos{k,t) + ^uj2ns{k,t), (36) 

n=l 

= / Qs(fc,t)dt + i V / r2„3(fc,t)gs(fc,0dt, (37) 

•^^s ^ „=i ^r^3 

- / Qs(fc,t)di + i V / /2„s(fc,i)di, 



(38) 



ujT{k,t) = ujQ^[k,t) + y^^2«T(fc:^): (39) 

n=l 

= [ QT{k,t)dt+l-y2 f Y2nAk,t)QT{k,t)dt, (40) 



^ QT(fc,i)dz + i^^ /2„^(fc,t)dt, 



where 



/2s(fc,t) 


= A{k,t){t- 




-5/2 


/4s(fc,i) 


= B{k,t){t- 


TsT 


-11/2 


h-rik.t) 


= C{k,t){t- 


ttY 


-5/2 


fi-vik^t) 


= D{k,t){t~ 


tt) 


-11/2 



(41) 



/2„s(fc,<) = r2„s(^.0Qs(A:,i), (42) 
f2nAk,t) - r2„^(fc,t)QT(fc,i). (43) 

The functions f2nsik,t) and f2nT{k,t) have the following functional dependence: 

(44) 
(45) 
(46) 
(47) 

where the functions A{k,t) and B(k,t) are regular at ts and the functions C{k,t), D{k,t) are regular at tt- With 
the help of the functions (|l^ - ([Tf|) we compute the integrals for uj2n up to = 4 using the contour indicated in 
Figs. |5(b)| -(c) and 6(b) (c). The expressions for u;2n permit one to obtain the fifth-order phase integral approximation 



of the solution to the equations for scalar (|16p and tensor (jl7p perturbations. The constants ci, C2, di and ^2 are 
obtained using the limit A; t ^ of the solutions on the left side of the turning point and (pij) , and are given by 
the expressions 



Cl = 


-1C2, 




e-'T 


C2 = 


— ^6 


di = 


-1^2, 




e-'T 


d2 = 


— 

V2 



-i[fcr,(0)-l-|(^Os(fc,0)|] 
-i[fej7(0)-f|(.c;oT(fc,0)|] 



(48) 

(49) 
(50) 
(51) 



In order to compute the scalar and tensor power spectra, we need to calculate the limit as /c i oo of the growing 
part of the solutions on the right side of the turning point given by Eq. ([32|) and Eq. ([33]) for scalar and tensor 
perturbations respectively. 



Ps{k) 
PAk) 



fc3 

hm —r 

-kt^OQ ZTT 



lim —r 



ait) 



(52) 
(53) 
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E. Uniform approximation 



We want to obtain an approximate solution to the differential equations (|16p and (fT7|) in the range where Qg(fc,t) 
and Q^{k,t) have a simple root at t^ct = t'Sj and <rot = ^t, respectively, so that Qgxl^?^) > for < t < irot 
and t(^; ^) < ^r t > <ret as depicted in Fig. 5(a) and Fig. 6(a) Using the uniform approximation method 
[H [Mill i^l , we obtain that for < t < i^t we have 



Uk{k,t) = 
Vk{k,t) = 



tret we have 

MM)_''^' 
QUKt) 



{CiM-pi{k,t)]+C2B,[-pi{k,t)]}, 



1/4 



{CiM-pi{k, t)] + C2B,[-py{k, t)]} , 

[Ql^{k,t)f^dt, 

" Jt 

where Ci and C2 are two constants to be determined with the help of the boundary conditions pTjl . For t > t^. 



(54) 
(55) 
(56) 



C/fe(fc,t) 
Ufe(fc,t) 



-Pr{k,t) 



-,1/4 



-Pr{k,t) 



{CiA,[p,{k,t)] + C2B,[p,{k,t)]}, 



1/4 



{CiA[p,{k,t)] + C2B,[p,{k,t)]}, 



13/2 _ 



r [-QiT{k,t)Y^^dt, 

Jtrnt 



(57) 
(58) 
(59) 



■J Jt 

For the computation of the power spectrum we need to take the limit kt ^ 00 of the solutions ([57| and (|58p . In 
this limit we have 



^hgi(M)]--^{iexp(-£[-gi(M)]-^d. 
iexp (^^* [-gi(fc,i)]'/'dt)| 

[-Q^(fc, t)] -^/^ I i exp £ [-Q|(fc, i)] dt 
iexp^* [~Q\{k,t)f\\i)jY 



(60) 



(61) 



where C is a phase factor. Notice that Eq. (pOj) and Eq. (pTj) are identical to Eq. ([5^ and Eq. obtained in the 
first-order phase-integral approximation. Using Eqs. pT|) . (fT3)) and the growing part of the solutions ([60)1 and ((6T|) 
one can compute the scalar and tensor power spectrum using the uniform approximation method, 



Ps(fc) = Hm ^ 



PT(fc) 



^3 

hm 



ait) 



(62) 
(63) 



We also use the second-order improved uniform approximation for the power spectrum [3], 



Ps,T(fc) = Ps,T(fc) [r*(i^S,T)], 

where I's.T is the turning point for the scalar or tensor power spectrum and 



(64) 



12t/ 288i^ 51840J/ 



(65) 



12 



F. Slow- roll approximation 



The scalar and tensor power spectra in the slow-roU approximation to second-order are given by the expressions 



Psik) 



1 + (4c ~ 2)ei + 2cA + ( 3c2 + 2c - 22 + e^Si + Uc^ - A +^] Sl- 



im) 



2^ 



k=aH 



1 + (2c - 2)ei + ( 2c^ - 2c - 3 + y 1 ef + ( -c^ + 2c - 2 + — ) £2 



where b is the Euler constant, 2 - hr 2 - 6 ~ 0.7296 and hr 2 + 6 - 1 ~ 0.2704, and 



(67) 



k=aH 



HI 

lp = TTm:' 

^1^ = ^ (69) 

' H dt (1 + 27V,)'' 
1 d"+^(l) 



The spectral index in the slow-roll approximation are 



n|'(fc) ~ l-4ei-2(5i-|-(8c-8)e? + (10c-6)eiJi, (71) 
n^'(fc) ~ -2ei-2e2 + (2a-2)e2. (72) 

The expressions (p7| . (1711) y (|7^ depend explicitly on time. In order to compute the scalar and tensor power 
spectra we need to obtain the dependence on the variable k. For a given value of k (0.0001 Mpc~^ <k< 15Mpc~^) 
we obtain t^, from the relation k = aH . Thus, for each k one obtains a value of t that we substitute into and Eqs. 

(IMl), (EZD, dZll) y 



G. Numerical solution 



We integrate on the physical time t the equations and (fT7|) governing the scalar and tensor perturbations using 
the predictor-corrector Adams method of order 12 [2^], and solve two differential equations, one for the real part 
and another for the imaginary part Uk and Vk- Two initial conditions are needed in each case UkiU), Uj^iU), VkiU), 
V^{ti), which can be obtained from the third-order phase- integral approximation. We start the numerical integration 
at tj calcu lated at 25 oscillations before reaching the turning point tret [SH. We call this procedure ICs phi3. Figs. 
7(a)||9(b) compare the numerical solution with the fifth-order phase-integral approximation for Re(ufc), lm{uk), \uk\, 



Ke{vk), lm{vk), and \vk\- Figures are plotted against the number of c-folds N. The solid line corresponds to the 
numerical solution (ICs phi3), the dashed line corresponds to the fifth-order phase- integral approximation. In each 
case the turning point rg and tt are indicated with an arrow. We stop the numerical computation of Psik) and Pt(^) 
a,t t = 5.00 X 10^ -^^pi^i after the mode leaves the horizon, where Uk/zg and Vk/a are approximately constant. Notice 
that the expressions for fitting pT|) and are valid in the aforementioned time scales, therefore we can use them 
for computing the scalar Psik) and tensor PT(fc) power spectra. 



H. Results 



For the chaotic ^m?'(fp' inflationary model, we want to compare the scalar and tensor power spectra and the 
spectral indices for different values of k calculated using the third and fifth-order phase-integral approximation with 
the numerical result (ICs phi3), the first and second-order slow-roll approximation and the first and second-order 
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FIG. 7: (a) Ke(uk) and (b) Re(i;fe) versus the number of e-folds for the chaotic inflationary model ^m^^^. Solid line: numerical 
solution (ICs phi3); dashed line : fifth-order phase-integral approximation 



uniform approximation method. First we analyze the resuhs for the scalar Ps{k) and tensor Prik) power spectra 
shown in Fig. [TO] and Fig. [T^ 

Table U shows the value of Ps{k), Prik), ns{k), and riT(fc) using each method of approximation at the WMAP 
pivot scale. It can be observed that the best value is obtained with the fifth-order phase-integral approximation. It 
should be noticed that the slow-roll approximation works well since the parameters ei, £2 and (5„ are small. 

Figures [TT] and [T51 show the relative error with respect to the numerical result that is obtained using the expression 
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FIG. 8: (a) Im(ttfc) and (b) Im(vfe) versus the number of e-folds for the chaotic inflationary model |m^</)^. Solid line : numerical 
result (ICs phi3); dashed line: fifth-order phase-integral approximation 



error rel. Ps,T(fc) = 



papprox 
-^S,T 



(fc) - Pg""™ 
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pnum 
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X 100. 



(73) 



The first-order phase integral approximation 
result, and deviate from the numerical result in 



the WKB and the first-order uniform approximation give the same 
10%. The second-order improved uniform approximation gives an 
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FIG. 9: (a) \uk\ and (b) \vk\ versus the number of de e-folds for the chaotic inflationary model 
result (ICs phi3); dashed line: flfth-order phase-integral approximation 



Solid line: Numerical 



error of 0.6%. With the first and second-order slow-roll approximation we have an error of 1% for Ps{k) and of 
0.4% for Prik). Using the third-order phase- integral approximation the error gives 0.15%, whereas the fifth-order 
phase-integral reduces to 0.015% in both cases. Fig. [Ml and Fig. [T5l show the results for the spectral indices ns{k) 
and nT(fc) respectively. 
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0.0001 0.001 0.01 0.1 1 10 

k[Mpc^] 

FIG. 10: Ps{k) for the chaotic inflationary ^m?(lP model. Thin sohd hne: numerical result (ICs phi3); dot-dashed line: third- 
order phase-integral approximation, dashed line: fifth-order phase-integral approximation: thick solid line: first-order phase- 
integral approximation, WKB and first-order uniform approximation, dashed double-dots line: second-order improved uniform 
approximation, double-dashed dot line: second-order slow-roll approximation, dotted line: first-order slow-roll approximation. 
The inset is an enlargement of the figure. 
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FIG. 11: Relative error for Ps(fc) for the chaotic inflationary ^m^cfy^ model. Dot-dashed line: third-order phase-integral 
approximation; dashed line: flfth-order phase-integral approximation; thick solid line: flrst-order phase-integral approximation; 
WKB and flrst-order uniform approximation; dashed double-dots line: second-order improved uniform approximation; dotted 
line: first and second-order slow-roll approximation. The inset is an enlargement of the figure. 
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FIG. 12: Prik) for the chaotic inflationary ^m?4P model. Thin solid lino: numerical result (ICs phi3); dot-dashed line: third- 
order phase-integral approximation; dashed line: fifth-order phase-integral approximation; thick solid line: first-order phase- 
integral approximation, WKB and first-order uniform approximation; dashed double-dots line: second-order improved uniform 
approximation; double-dashed dot line: second-order slow-roll approximation; dotted line: first-order slow-roll approximation. 
The inset is an enlargement of the figure. 
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FIG. 13: Relative error for ft'(fc) for the chaotic inflationary im^</>^ model. Dot-dashed line: third-order phase-integral 
approximation, dashed line: fifth-order phase-integral approximation; thick solid line: first-order phase-integral approximation; 
WKB and first-order uniform approximation; dashed double-dots line: second-order improved uniform approximation; dotted 
line: first and second-order slow-roll approximation. The inset is an enlargement of the figure. 
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FIG. 14: ns(fc) for the chaotic inflationary ^m'^(f>'^ model. Thin solid line: numerical result (ICs phi3) and fifth-order phase- 
integral approximation; dot-dashed line: third-order phase-integral approxiniation;tliick solid line: first-order phase-integral 
approximation; WKB and first-order uniform approximation; dashed double-dots line: second-order improved uniform approx- 
imation; double-dashed dot line: second-order slow-roll approximation; dotted line: first-order slow-roll approximation. The 
inset is an enlargement of the figure. 
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FIG. 15: nT(fc) for the chaotic inflationary ^m'^(f>'^ model. Thin solid line: numerical result (ICs phi3) and fifth-order phase- 
integral approximation; dot-dashed line: third-order phase-integral approximation; thick solid line: first-order phase-integral 
approximation; WKB and first-order uniform approximation; dashed double-dots line: second-order improved uniform approx- 
imation; double-dashed dot line: second-order slow-roll approximation; dotted line: first-order slow-roll approximation. The 
inset is an enlargement of the figure. 



19 





num" 


phi3'' 




srl" 


sr2'= 


phil'', WKB«, ual" 


ua2' 


Ps(fc) X 10-" 


8.3300 


8.3419 


8.3287 


8.2282 


8.2299 


7.4946 


8.2752 


Pt(A:) X 10-1^ 


1.6109 


1.6132 


1.6106 


1.6047 


1.6049 


1.4486 


1.6004 


ns(fc) 


0.960473 


0.960472 


0.96473 


0.960700 


0.960491 


0.960453 


0.960486 


mik) 


-0.019965 


-0.019982 


-0.019965 


-0.019650 


-0.019948 


-0.020000 


-0.019984 



TABLE I: Value of Ps(k), PT{k), ns{k), and nT(fc) obtained with different approximation methods for the chaotic inflationary 
model ^m'^ff'^ for the mode k = 0.05 Mpc~^. 

"Numerical 

''Third-order phase-integral approximation 
"Fifth-order phase-integral approximation 
■"First-order slow-roll approximation 
"Second-order slow-roll approximation 
'First-order phase integral approximation 
*WKB approximation 
'Tirst-order uniform approximation 
'Second-order improved uniform approximation 



III. CONCLUDING REMARKS 



The results reported in this article show that, in comparison with other approximation methods, the phase integral 
approach gives very good resuhs for the scalar and tensor spectra in the quadratic inflationary model. The phase- 
integral approximation gives very accurate results as soon as the the integral /^(z, zq) is small. Figures 7(a)||9(b") 



show that the phase integral approximation fails in the vicinity of the turning point —v, range where the /i-integral 
diverges. The selection of the base function Q{z) guarantees that /i <C 1 far from the turning point at any order of 
approximation. Since the scalar and tensor power spectra as well as the spectral indices are evaluated as —krj—f 0, 
the limit is taken far from the horizon (turning point), therefore their computation is not affected by the presence of 
the turning point. 

Since the WKB method can be regarded as a first-order approximation of the phase-integral approximation with 
Q^{z) = R{z), it should be expected that the phase- integral method works in those cases where the WKB methods 
give s g ood estimates and slow-roll fails, that is the case where inflation is generated by a chaotic potential with a step 
[ISl . |33 | . The good agreement between the numerical results and those obtained with the phase-integral approximation 
shows that the phase integral method is a very useful approximation tool for computing the scalar and tensor the 
power spectra in a wide range of inflationary scenarios. 
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